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Using data from both solar and atmospheric neutrinos, Super-Kamiokande has provided fun- 
damental information on neutrino flavor mixing and neutrino mass square differences. 
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1 Introduction 

X 

Both solar and atmospheric neutrino measurements have yielded significantly smaller interaction 
rates than expected: the reduction of tl^ solar u e flux compared to the standard solar model 
(SSMn) was observed to be between 50%BBl3ti and 65% 00, the atmospheric v e flux was found to 
be the same within uncertainty while the atmospheric flux was measured at about 65% of the 
calculation^. Both deficits — known as the solar neutrino problem and the atmospheric anomaly 
- are explained by neutrino flavor changes caused by oscillations: solar neutrino experiments 
detect predominantly v e and atmospheric experiments are less sensitive to v T interactions. 

Neutrino oscillations require the physical mass eigenstates z/j (with masses rrt\ < ni2 <■■■) 
to differ from the flavor eigenstates v a (defined by charged-current interactions with a charged 
lepton of flavor a). The flavor eigenstates v a are orthogonal linear combinations of the mass 
eigenstates f^: v a = ^ U a iVi. U is the unitary mixing matrix. Since there are three charged 
lepton flavors (e, [i and r), only three neutrino flavor eigenstates can produce charged leptons 
via charged-current interactions. The invisible width of the Z° boson is compatible with decay 
into three light neutrinos, so only three neutrino flavor eigenstates participate in neutral-current 
interactions. If there are more than three neutrinos (and therefore more than three neutrino 
flavor eigenstates), then some of the flavor eigenstates are "sterile", that is, they feel neither 
charged- nor neutral-current interactions. In case of just three neutrinos (and neglecting CP 
violating phases), the mixing matrix U can be described phenomenologically using the three 
mixing angles #12 , #13 and 623 ■ Using the short notation = sin 9{j , Cij = cos 6ij the mixing 



matrix becomes 
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and the flavor conversion probability for a neutrino of energy E after flight distance L in vacuum 
is 
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Atmospheric neutrino data requires! a Am^ m around 10 3 eV 2 , while solar neutrino data 
implies a Am^j below rs 10~ 4 eV 2 <C Am^ m . For the purpose of this paper, the Am^j 
is assumed to be Amf 2 - Then, Am^ m Amj 3 



Am 2 ^ and the atmospheric conversion 



probabilities are 
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while the solar i/ e survival probability is 
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with 6* so j = 6*12 and sin0 a t m = cos 6*13 sin 6*23 . Therefore, atmospheric neutrinos yield informa- 
tion about the last column of the mixing matrix and solar neutrinos about the first row. If 6*13 
is small, then both solar and atmospheric neutrinos can be analyzed using a two-neutrino model 
with effective solar and atmospheric mixing angles. 

The characteristic L/E dependence of the vacuum oscillation probability is modified by the 
matter density along the flight path of the neutrino, since the coherent scattering amplitude 
of v e with matter is different from Vn or v T , and sterile neutrinos don't undergo such coherent 
scattering at all. To take into account these matter effects, the potential V can be defined. The 
matter potential is the same for j/„ and u T , but differs for v e (or a sterile flavor). The same 
two-neutrino formulas may be used, with the substitution of 
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for Am 2 and sin 2 29, where AV is the potential difference. While the vacuum conversion prob- 
abilities are invariant under 9 — ► n/2 — 9 (that is sin 20 — ► sin 20) the cos term changes sign 
(cos 26* — ► — cos 29) and the effective Am 2 and mixing changes. Therefore the "mirror symme- 
try" around maximal mixing (9 = vr/4) is broken by the matter effects. Depending on the sign 
of AV, a resonance can occura which leads to a large conversion probability, even if the vacuum 
mixing is small. 
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Figure 1: SK Atmospheric Neutrino Data and Atmospheric Neutrino Oscillation Parameters, (a) The upper 
four panels show the zenith angle distribution of single ring events for fully contained e-like rings (panel one and 
three), fully contained /i-like rings (panel two) and higher energy /i-like rings (panel four). In the lower row, 
panel one (brightest ring is /i-like) and two (brightest ring is e-like: neutral-current enhanced sample) contain the 
zenith angle distributions of the multi-ring samples. The zenith angle distributions of interacting in the rock 
below SK are given in panel three (lower energy) and four (higher energy). The predictions of the atmospheric 
neutrino Monte Carlo calculation are overlaid (light-gray hatched areas). The best oscillation fit is indicated 
by the dark-gray solid line, (b) Allowed oscillation parameters for a two-neutrino fit assuming pure i/ M <-» v T 

oscillations. 



2 Atmospheric Neutrinos 

Atmospheric neutrinos originate when primary cosmic rays strike the earth's atmosphere and 
produce a shower of pions, which decay into muons and z^'s. Each muon decay produces a u e 
as well as a v^. Super-Kamiokande (SK) is a cylindrical 50,000 ton water Cherenkov detector. 
An optical barrier separates an "inner" concentric cylinder from the outer "anti-counter" . The 
inner cylinder (32,000 tons) is viewed by 11,146 inward-facing 20" photomultiplier tubes, the 
outer cylinder by 1885 outward-facing 8" photomultiplier tubes. SK uses the anti-counter (which 
surrounds the inner detector) to define "fully contained events" (no entering or exiting charged 
particles), "partially contained events" (only exiting charged particles) and "entering events". 
In each case, the direction of a charged particle is reconstructed using the directionality of the 
Cherenkov light. 

Cherenkov rings are predominantly produced by electrons and muons. The Cherenkov ring 
produced by electrons typically has a "fuzzy" outer edge due to electromagnetic showers. A 
muon, on the other hand, does not shower as much, and therefore its Cherenkov ring has a 
"sharp" outer edge. This allows the construction of a particle identification likelihood (PID) 
which distinguishes "e-like" from "/i-like" events. The energy of fully contained events is mea- 
sured using the amount of produced Cherenkov light. The sample is split into a "sub-GeV" and 
a "multi-GeV" sample. If the brightest ring in an event containing several Cherenkov rings is 
/x-like, the event is most likely due to a charged-current interaction. Since the decay photons 
of neutral pions produce electrons by Compton scattering, a neutral-current enhanced multi-ring 
sample can be defined using the PID and the 7r° invariant mass. 

Partially contained events are almost always due to higher energy muons. The energy of 
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Figure 2: Neutrino Mixing with Three Parameters, (a) Allowed area of sin 2 #13 and sin 2 823. At 90% C.L. SK data 
allows up to about 40% of e-type content of the atmospheric oscillation, (b) Allowed range of the atmospheric 
v e disappearance probability compared to the limits for Z^T obtained from various reactor experiments. The most 

stringent limit comes from the CHOOZ experiment. 



those muons can only be estimated statistically using Monte Carlo distributions. Upward-going 
entering events are due to the interaction of i/„ in the rock below SK. If the "upward-going 
muons" stop inside the detector ("upward/stopped"), then the average energy of the producing 
v^'s is lower than for upward-going muons which enter and exit the detector ("upward/through"). 

The direction of low energy leptons yield little information about the neutrino's direction 
(and in turn its flight distance). However, the zenith angle 6 of high energy leptons is correlated 
with the neutrino flight distance L. SK has accumulated 1289 live days of atmospheric neutrino 
data using a fiducial volume of 22,500 tons. The zenith distributions of electrons (see figure [[]) 
show good agreement with the Monte Carlo based on the atmospheric neutrino flux calculation. 
However, there are significantly less muons than expected. Downward-going (cos 8 = 1) atmo- 
spheric neutrino-induced (multi-GeV) //-like events agree with the Monte Carlo, upward-going 
events show the strongest deficit. A similar distortion is seen in the multi-ring /i-like sample, 
but not in the multi-ring neutral-current enhanced sample. Upward-going through-going muons 
agree with Monte Carlo in the horizontal direction, but are suppressed at cos 6 = — 1. There are 
less upward-going stopped muons than expected. As seen in figure [j] (a) oscillations of u^s into 
z; r 's are able to explain all the zenith angle distributions. The same figure also gives allowed 
ranges for neutrino mixing and Am 2 . The data prefer maximal mixing and Am 2 ~2-3-10~ 3 eV 2 . 

Due to the lack of an e-like appearance signal, <-> u e oscillations are disfavored compared 
to <-> v T . However, any e-like appearance due to Vn —* v e is unfortunately washed out by 
a disappearance due to v e — > z/„, so the SK data allow a substantial V e " component (about 
40% at 90% C.L.) in the oscillations. Figure |2| shows the allowed values for both #13 and #23 
(according to equation (jl]) Qyi plays no role for atmospheric neutrino oscillations) in a three- 
neutrino oscillation analysis. Maximal mixing for #23 and zero mixing for #13 is preferred. The 
same figure also gives allowed ranges for Am^ m and the "disappearance probability" sin 2 2^13. 
More stringent limits on this probability come from the CHOOZ reactor experiment 0. 

Four-neutrino oscillations are investigated for the case of #13 = and a Am 2 hierarchy 
(Am<L SND « leV 2 > Am 2 tm « 10~ 3 eV 2 > Ara^j « 10~ 4 eV 2 ). The "content" of the 
fourth flavor (sterile neutrino) in atmospheric neutrino oscillations can be parameterized by 




Figure 3: Limit on Sterile Content of Atmospheric Neutrino Oscillations, (a) Pure <-> v T (solid black line) fits 
the upward-going fi zenith distribution better than pure z/ M <-» v s (solid gray line) because of matter effects. For 
the same reason, there is also a small difference in the high energy /^-like events. The neutral-current enhanced 
multi-ring sample shows no hint for a disappearance of active flavors and therefore favors pure <r+v T , (b) The 

limit on the sterile content is sin 2 £ < 0.25 at 90% C.L. 



sin 2 £: — > (cos£v T + sm^u s ). Figure ||| compares the best pure — > v T fit with the best 
pure —s- v s fit and gives the allowed ranges for Aro^ and sin 2 £. In case of purely sterile 
oscillations, matter effects reduce the conversion probability for high energy neutrinos which 
affects the upward-going muon samples, the partially contained events, and the multi-GeV [i- 
like events. The neutral-current enhanced multi-ring zenith angle distribution becomes up/down 
asymmetric, since sterile neutrinos don't interact via neutral currents. In every case, data 
disfavors the purely sterile case. Data likes best zero sterile content At 90% C.L., the maximum 
allowed sterile content is 25%. The range of Am 2 is consistent with the other fits. 

Other hypotheses to explain the SK zenith angle distributions are investigated as well. 
The conversion probability is parameterized as follows: (i) sin 2 29 sin 2 (aL x E) ("LxE"), (ii) 

sin 4 9 + cos 4 (l - e" aL / £ ) decay") and (iii) (sin 2 9 + cos 2 9e~ aL / E ^ 2 (V M decay to v„"). 
Table |l| gives a summary of all investigated hypotheses, — > v T is strongly favored by the data. 

Table 1: Comparison of Several Flavor Conversion Hypotheses, — > v T is strongly favored by the data. 



Mode 


best fit parameters 


Ax 2 


#^ 


Up -> u T 


sin^ 29 


= 1.00 


Am' 2 = 2.5 • lO^eV 2 


0.0 


0.0 


->• v e 


sin 2 29 


= 0.97 


Am 2 = 5 • 10~ 3 eV 2 


79.3 


8.9 


v^i -> v s 


sin 2 29 


= 0.96 


Am 2 = 3.6 • 10" 3 eV 2 


19.0 


4.4 


LxE 


sin 2 29 


= 0.90 


a = 5.3 • 10- 4 /(GeVkm) 


67.1 


8.2 


decay 


cos 2 9 


= 0.47 


a = 3.0 • 10" 3 GeV/km 


81.1 


9.0 


i/ju decay to u s 


cos 2 9 


= 0.33 


a = 1.1 • 10" 2 GeV/km 


14.1 


3.8 



If atmospheric v„ indeed oscillate into v T , then about 80 v T charged-current events are 
expected in the SK data set. The high threshold for this reaction makes it difficult to isolate 
these events. At this energy, atmospheric neutrinos typically produce a spray of particles, 
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Figure 4: Search for Tau Appearance, (a) Neural net output for atmospheric charged-current interaction Monte 
Carlo (dotted line), all atmospheric v Monte Carlo without v T (dashed line) and tau Monte Carlo (solid line). 
Since there is only a small number of tau events expected in the SK data, the analysis is limited by background 
from other atmospheric //-induced events, (b) Fit to the zenith angle distribution of the tau-enriched sample. Tau 
events are only expected in the upward direction. Data favors tau appearance at about the 2<r level. 



and ring-counting becomes very difficult. One analysis abandons ring-counting altogether and 
reconstructs the amount of "energy flow" as a function of direction based on the angular light 
distribution. From the energy flow, jet variables are formed. A r likelihood function combines 
such variables to discriminate between r events and other atmospheric neutrino interactions. 
Another analysis constructs a similar r likelihood from variables based on the conventional 
event reconstruction (ring-counting) A third analysis utilizes a neural net to combine similar 
ring-counting variables. Figure ^| (a) explains the output of the "tau neuron" when tested on 
Monte Carlo. In spite of the good separation between tau and non-tau events, the analysis (like 
the other two) is limited by background, since there are much fewer tau events than non-tau 
events. The same figure shows the zenith angle distribution of the tau-enriched sample (using 
the tau neuron). An excess of about two sigma (see also table ||) is observed for upward-going 
events. 



Table 2: Results of Three r Appearance Searches. 



Analysis 


Number of r-Events in Fit Efficiency 


Observed, Expected Signif. 


Energy-Flow 
Ring-Counting 
Neural Net 


79+^(stat+syst) 32% 
66±41 (states (syst) 43% 
92±35(stat)+23(syst) 51% 


1.8cr 1.9a 
1.5a 2.0a 
2.2a 2.0a 



3 Solar Neutrinos 

SK observes solar neutrinos via elastic neutrino-electron scattering. Only solar neutrinos with 
fairly high energy (several MeV) can be detected. These neutrinos originate from the /3 + -decay 
of 8 B or the 3 He-proton (hep) fusion reaction. Like other experiments SK found a deficit of solar 
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Figure 5: Time Variation of the Neutrino Flux Inferred from the Elastic Scattering Rate. The black points show 
the elastic scattering rate in SK during 1.5 months. The gray points extend over a time period of 10 days. Each 
black point is significantly below the SSM prediction and above the SNO charged-current rate. The solid line 
shows the 7% flux variation which is expected from the 3.5% change in the distance between the sun and the 

earth. 



neutrinos; all such deficits are explained by solar neutrino oscillations. Since event times are 
recorded and the event energy is reconstructed, SK can search for time variations and an energy 
dependence of the rate deficit. 

The solar neutrino flux is inferred from the SK interaction rate assuming only v e J s. The 
large solar neutrino event sample collected by SK during 1496 live days inside a fiducial volume 
of 22,500 tons enables a precise study of the time-dependence of this flux. The gray data points 
in figure [| contain only 10 days of data. The flux measured in each of these 185 time bins 
is significantly smaller than the SSM prediction. Also shown are the larger 1.5 month time 
bins (black data points). At first glance, the flux looks constant. A more careful study reveals 
an indication for a periodic change with a one year period. Since the sun is a neutrino point 
source, a 7% change of that kind is expected by the 3.5% change in distance between the sun 
and the earth (see solid black line of figure [5]). If the 1.5 month data points are combined into 8 
"seasonal" bins, a x 2 test to this 7% variation yields x 2 = 4.7 (69% C.L.) compared to x 2 = 10.3 
(17% C.L.) for a constant flux. The difference in x 2 corresponds to about 2.5 sigma. There is 
no hint of other long-term time variations (e.g. arising from the oscillation phase). 

The charged-current interaction rate measured by the SNO experiment explores the solar 
neutrino flux in a similar energy range and can therefore be compared directly. The flux mea- 
sured by SK in each of the 42 (1.5 month) time bins is significantly larger then the flux measured 
by SNO. Solar neutrino oscillations of v e into other, active flavors predict such an excess of the 
elastic scattering rate over a charged-current measurement; that excess is due to a larger neutral- 
current amplitude in the elastic-scattering process. The discrepancy between the SK and the 
SNO neutrino flux measurement is therefore interpreted as evidence for the presence of other 
active flavors than v e in the solar neutrino flux. 




Figure 6: Zenith Spectrum and Four Oscillation Predictions. SMA (dotted line) and VAC (dashed line) solutions 
are disfavored by spectral distortion data, LOW (light gray) solutions by the absence of solar zenith angle variation. 

The LMA solutions (dark gray) fit best. 



Depending on the solar Am 2 , solar neutrino oscillations can be driven by the oscillation 
phase or by matter effects. Strong spectral distortions are expected in the former case. In 
the latter case, spectral distortions can still occur as a consequence of a (energy-dependent) 
resonance caused by the sun's matter density. The earth's matter density can also affect the 
conversion probability resulting in most cases in a "regeneration" of "disappeared" u e 's during 
the night. SK searches for both effects with the "zenith angle spectrum" (see figure ^): The 
data is broken into eight energy bins (according to the SK energy resolution). Between 5.5 and 
16 MeV of reconstructed energy, SK has collected a sufficiently large number of solar neutrino 
candidate events to subdivide each energy bins in seven solar zenith angle (6 Z ) bins. In each 
bin, the ratio of the observed event rate and the rate expected by the SSM is shown. Neither 
spectral distortion nor a regeneration effect is evident in the data. 

For the oscillation analysis ffl of this zenith angle spectrum (and other solar data), #13 is 
set to zero in equations @. This is justified by the lack of evidence for a positive $13 in either 
atmospheric neutrino data from SK or reactor neutrino data from CHOOZ. Analysis of the 
data from the other SK experiments reveal several distinct allowed areas shown in light-gray in 
figure |7| (a) which are referred to as "solutions" to the solar neutrino problem. The small mixing 
angle solution (SMA) is the only one with a small solar mixing angle. There are three large angle 
solutions. The large mixing angle solution (LMA) has the largest Am 2 and the vacuum solution 
(VAC) the smallest. The lower part of the region extended between 10 _9 eV 2 and 10 _7 eV 2 is 
called quasi- VAC solution, the upper part is the LOW solution. Figure ^ shows the typical 
spectral and solar zenith angle behavior of the main four solutions. The absence of spectral 
distortions excludes VAC and SMA, the absence of the nightly "regeneration shape" required 
by the LOW solution disfavors this solution as well. Figure [?] (a) shows the excluded regions of 
parameter space due to only the shape of the zenith angle spectrum (gray areas). Only LMA 




Figure 7: Solar Neutrino Oscillation Parameters, (a) The gray area is excluded at 95% C.L. by the shape of the 
SK zenith angle spectrum, the light-gray area is allowed by data from Gallex/GNO, SAGE, Homestake and SNO. 
The overlap of both is shaded dark-gray, (b) The light-gray area is allowed by the shape of the SK zenith angle 
spectrum and the SK interaction rate using the SSM prediction of the 8 B flux, (c) The light-gray area is allowed 
using all solar neutrino data. The SSM predictions of the 8 B and hep flux are not used. The best fit (located in 

the LMA solution) is marked by the white asterisk. 



and quasi- VAC solutions remain allowed. Combining the shape of the zenith angle spectrum 
with the SK rate measurement, only the two allowed regions in figure |7] (b) remain (LMA and 
quasi-VAC). SK solar neutrino data by itself requires large mixing and prefers maximal mixing. 
Stronger constrains are obtained when the data of the other solar neutrino experiments are used 
in addition to SK data. The quasi-VAC solution is rejected by the rates of all experiments. 
Combining £3 SK zenith angle spectrum and rate with the rates of all other solar neutrino 
experiments, SK found for the first time a unique solution, the LMA solution, at ~95% C.L. 
Figure ^ (c) shows the allowed range in mixing and Am 2 . The mixing is large but not maximal 
and Am 2 is somewhat smaller then 10 -4 eV 2 . Since SK spectral data limits high energy solar 
neutrino spectral distortions to be small, the combination of the SK rate and the charged-current 
rate of either the Homestake experiment or SNO forbids a purely sterile solar neutrino oscillation 
solution. 



4 Neutrino Mixing and Masses 



Figure || summarizes and combines the current best knowledge of neutrino masses and mixing 
from solar and atmospheric neutrino data. There are no hints for sterile neutrinos in either 
atmospheric or solar neutrino oscillations, so three neutrinos are sufficient to explain both data 
sets. The mass spectrum of figure |§| is for the hierarchical case, that is m% is much smaller than 
the other masses, and the mass eigenstate pair responsible for solar neutrino oscillations is m\ 
and m2. CP violating phases are neglected. The displayed values of Oij are the best-fit values 
from the atmospheric and the solar oscillation analyses. 



u 



' 2 

1_ 

1 

2V2 



712 « fi 



723 



1 

2 

2^ 
1_ 

2y/2 



\ 

1 

72 / 



4 #13 







$ 50 
E 

c 40 

S3 30 
E 



20 
10 



e-type content 
l^-type content 
x-type content 



m 



Figure 8: A Possible Leptonic Mixing Matrix U and Neutrino Mass Spectrum. The mixing matrix neglects 
CP violation and assumes that #13 is indeed zero and the atmospheric mixing is maximal. The mass spectrum 
additionally assumes, that mi is small enough to be negligible and that the masses of the "solar pair" are smaller 

than the third mass, (hierarchical case). 
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